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Coevolution Between Supermassive Black Holes and Bulges Is Not Via 
Internal Feedback Regulation But By Rationed Gas Supply Due To 

Angular Momentum Distribution 

Renyue Cen^ 

ABSTRACT 

We reason that, without physical fine-tuning, neither the supermassive black 
holes (SMBHs) nor the stellar bulges can self-regulate or inter-regulate by driving 
away already fallen cold gas to produce the observed correlation between them. 

We suggest an alternative scenario where the observed mass ratios of the SMBHs 
to bulges reflect the angular momentum distribution of infallen gas such that the 
mass reaching the stable accretion disc is a small fraction of that reaching the 
bulge region, averaged over the cosmological time scales. We test this scenario 
using high resolution, large-scale cosmological hydrodynamic simulations (without 
AGN feedback), assuming the angular momentum distribution of gas landing in the 
bulge region to yield a Mestel disc that is supported by independent simulations 
resolving the Bondi radii of SMBHs. A mass ratio of 0.1 - 0.3% between the very 
low angular momentum gas that free-falls to the sub-parsec region to accrete to the 
SMBH and the overall star formation rate is found. This ratio is found to increase 
with increasing redshift to within a factor of ~ 2, suggesting that the SMBH to bulge 
ratio is nearly redshift independent, with a modest increase with redshift, a testable 
prediction. Furthermore, the duty cycle of active galactic nuclei (AGN) with high 
Eddington ratios is expected to increase significantly with redshift. Finally, while 
SMBHs and bulges are found to coevolve on ~ 30 - ISOMyr time scales or longer, 
there is indication that, on shorer time scales, the SMBH accretion rate and star 
formation may be less correlated. 


1. Introduction 


There is mounting evidence that massive bulges in the nearby universe harbor central 
SMBHs of mass 10® - 10® M©. The correlation between SMBH mass (Mbh) and the bulge 


(BG) mass ( Mbg) or velocity dispersion (a) (e.g., Magorrian et al. 

1998;[Richstone et al. 

1998 

Gebhardt et al. 2000;[Ferrarese & Merritt 2000 Tremaine et al. 

20021 suggests coevolution. 


Although alternative models for producing this observed relation are available (e.g., [Ostriker 
20001 [Adams et al.|2001t[Colgate et al.||2003t|Cen|2007| , the correlation is often construed as 
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evidence for AGN feedback to regulate the growth of SMBHs and bulges. The idea that AGN 


feedback may alleviate problems in galaxy formation models (e.g., Kauffmann & Haehnelt 


20001 [Croton et'ar]|2006t jSomerville et al.|[200^ further enhances its appeal. The three- 
dimensional hydrodynamic simulations successfully reproduced the observed Mbh/Mbg ratio 
(e.g., |Di Matteo et al. 2005[|Hopkins et al.|2006| , providing the physical basis for this scenario. 


This Letter has two goals. First, we make a qualitative examination of the implications of 
the observed relation between bulges and the central massive objects (CMOs), wherein the 
two follow a linear relation over four decades in mass. It is shown that neither the SMBHs nor 
the nuclear star clusters (NSCs) nor the stellar bulges could have played a dominant role in 
regulating the growth of any of the three components in the way of blowing away a significant 
fraction of gas already landed in the respective regions so as to produce the CMO-bulge 
relation. Second, an alternative model is put forth wherein the correlation between SMBH 
mass and bulge mass is dictated by the angular momentum distribution of the infalling gas. 
We successfully test this new scenario using ab initio Large-scale Adaptive-mesh-refinement 
Omniscient Zoom-In (LAOZI) cosmological hydrodynamic simulations. 


2. Arguments Against Internal Regulation of the Central Components 


With the ACS Virgo Cluster Survey of early-type galaxies spanning four decades in mass, 
Cote et alTj ( |2006| and [Ferrarese et al.] ( [2006^ find a transition at Mb,o = -20.5, where the 
brighter galaxies lack resolved stellar nuclei and SMBHs dominate the CMO mass, while fainter 
ones have resolved stellar nuclei that dominate the CMO mass. Furthermore, the logarithm of 
the mean nucleus-to-galaxy luminosity ratio in fainter, nucleated galaxies, -2.49 ± -0.09 (a = 
0.59 ± -0.10) is indistinguishable from that of the SMBH-to-bulge mass ratio, -2.61 ± -0.07 
{a = 0.45 ± -0.09). A similar result is found by |Wehner & Harris| ( |2006^ using a different data 
set. [Turner et ar||2012| find an identical relation using early-type galaxies in the ACS Fornax 
Cluster Survey. We express the universal scaling relation between CMOs and bulges as 


Mcmo — Mbh + Mnsc — Oi Mbg, 


( 1 ) 


whereby with the transition between NSC and SMBH occurs at Mb ~ —20.5 or stellar mass 
Mbgo = (3~4) X10^° Mq, and a = 2.5 x 10“^ |C6te et ar2006 ; Ferrarese et al. 2006 : Wehner 
& Harris[[2006t[Turner et al.[[2012| . One may express regulation of the growth of bulges as 


gbh Mbh + gnsc Mnsg + gbg Mbg = Mbg, 


( 2 ) 


where cbr, gnsc and bbg are the feedback strength coefficients per unit mass of the respective 
components exerted on the stellar bulge and the ejected gas mass is equal to IMbg; cr is the 
velocity dispersion of the stellar bulge; /3 is a parameter that absorbs uncertainties regarding 
the dynamics of concerned feedback processes, with /3 = 2 for energy-conserving feedback 
(gbh, eNsc and cbg have units of energy per unit mass) and /3 = 1 for momentum-conserving 
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feedback (eeH, gnsc and eec have units of momentum per unit mass). Note that a significant 
feedback regulation means / > 1. 

Insights can be gained by asking the following question: Can the feedback from SMBH 
and NSCs conspire to regulate the growth of the stellar bulge, i.e., 

eBH Mbh + gnsc Mnsc = fcr^ Mbg? (3) 


The single powerlaw relation between Mcmo and Mbg across four decades in bulge mass 
can be understood, only if the negative feedback per unit stellar mass of the NSC and of the 
SMBH are approximately the same, cbm ~ gnsg, barring the unknown physical reason for the 
right hand side of Eq ([^ - the required amount of notional feedback to regulate the bulge 
growth - to change character abruptly at Mbg — Mbgo- 

Although having bbe ~ gnsg niay be possible, it would render a negative answer to the 
question above (Eq[^, as follows. In the momentum driven regime, since the feedback from 
the nuclear cluster is subject to higher densities and shorter cooling timescales hence dimin¬ 
ished strength in comparison to that in the stellar bulge, i.e., bbg > eNsc- In the energy driven 
feedback scenario, bbg = gnsc- Since Mnsc < Mbg, the supernova feedback from stars in 
the bulge would vastly exceed that from the NSC. This thus invalidates the statement that the 
NSC and SMBH provide the necessary feedback to regulate the growth of the bulge. 


The only scenario left for the SMBH to regulate the bulge growth is to force cnsc = 0 
and assume the feedback per unit SMBH mass, while constant at Mbg > Mbgo, to become 
negligible at about Mbg = Mbgo- In both the momentum (/3 = 1, jOstriker et al.||2010| ) 
and energy feedback scenario (/? = 2, [Faucher-Giguere & Quataert||2012| ), the amount of 
momentum or energy per unit SMBH mass, cbb, is ultimately proportional to the driving energy 
(oc Mbhc^, where c is speed of light). Thus, there exists no known process to suddenly make 
eBH drop to zero at some specific Mbh, while being constant otherwise. 


If negative feedback is needed to internally regulate the bulge, the only alternative left is 
stellar feedback from bulge stars themselves, i.e.. 


Bbg = 


(4) 


Under the assumption that the feedback strength from stars per unit mass (cbc) is constant, 
one obtains / oc which has the same dependence on a as the predicted mass load¬ 
ing factors for both momentum (/3 = 1) or energy (/3 = 2) driven winds (e.g., [Murray et’aT 
2005| ). Therefore, bulge self-regulation, if required, would be physically supportable and self- 
consistent. If bulge is self-regulated, then, under the assumption that eNsc = gbg, NSC may 
also be self-regulated. The correlation between Mmco and Mbg would then require that the 
mass loading factor for the SMBH is the same as for the NSC, i.e., bbh = gnsc, which is a 
fine-tuned outcome. In the absence of inter-regulation between CMOs and bulges, the propor¬ 
tions of the amount of gas feeding the nuclear and bulge regions must be proportional to the 
observed Mcmo/Mbg I'stio. 
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3. An Alternative Scenario: Rationed Cold Gas Supply to Nuclear and Bulge Regions 

Over Cosmological Time Scales 

Our arguments in the previous section indicate that the observed Mcmo-Mbg correlation 
requires the same proportionality in the initial amounts of gas feeding the respective regions, 
averaged over the cosmological time scales. We test this scenario using direct cosmological 
simulations. 


3.1. Simulation Characteristics 


See |Cen||2014| for a more detailed description of the ab initio LAOZI simulations. Briefly, 
we use the WMAP7-normalized ( [Komatsu et al.||201 1 1 ) ACDM model: Qm = 0.28, Qb = 0.046, 
= 0.72, as = 0.82, Hq = 100/ikms“^Mpc“^ = 70kms“^Mpc“^ and n = 0.96. A zoom-in 
box of size 21 x 24 x 20/i“^Mpc^ comoving is embedded in a 120 /i“^Mpc periodic box. The 
maximum resolution is better than 111/i'^pc (physical) at all times. Star formation follows the 
prescription of jCen & Ostriker| ( [T992| . Supernova feedback from star formation is modeled fol- 
lowing [Ceh et al.| ( [20^ with feedback energy being distributed into 27 local gas cells weighted 
by the specific volume of each cell, to mimic the process of supernova blastwave propagation 
to channel more energy into the less dense regions. We exclude AGN feedback in order to 
ascertain the lack of need for it. 


3.2. Construction of Gas Feeding Histories of Simulated Galaxies 


Galaxies are identified using the HOP algorithm ( jEisenstein & Hut||1998[ ) grouping stellar 
particles. Galaxy catalogs are constructed from z = 0.62 to z = 1.40 with an increment of 
Az = 0.02 and from = 1.40 to x = 6 with Ax = 0.05, having a temporal resolution of 
30 - 150Myr. For each galaxy at x = 0.62 a genealogical line is constructed up to x = 6, where 
the parent of each galaxy is identified with the one at the next higher redshift with the most 
overlap in stellar mass. At each redshift, we compute the amount (Me) and mean specific 
angular momentum (Jc) of gas in the central Ikpc region. To proceed, an ansatz is made: 
the gas mass with angular momentum lower than Jn is Mc(/3Jn/(l + /9)Jc)^- We use (3 = 1 , 
which corresponds to a|Mestel|(|l963| disc of surface density S(r) oc (3 = l\s motivated 


by simulations of [Hopkins & Quata^ ( [2010[|2011) with resolution as high as 0.1 pc. Figure 
12 of [Hopkins & Quata^ ( [2010^ shows that the evolved density runs of the gas discs, on 
average, follow the S(r) oc r~^ profile from 0.1 pc to 1 kpc. In all of the six individual cases with 
significant gas inflow, shown in Figures (2, 3) of [Hopkins & Quatae7t|(|201 1|, the S(r) oc r~^ 


profile provides an excellent fit. We compute the 1-d stellar velocity dispersion a within the 
effective radius for each galaxy in the simulation at any redshift and assume an SMBH of 
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mass equal to Mbh = lO®M0((j/2OOkm/s)'‘ ( [Tremaine et al.||2002[ ). The Bondi radius is 

tb = 2GMBH/3cr^ = 7.2pc((T/200 km/s)^, (5) 

and the specific angular momentum at ib is 

Jb = V2TBa. (6) 


The gas landing within iq is assumed to accrete to the SMBH, where at r > iq the disc has 
Toomre Q parameter below unity and is hence consumed by star formation. Expressing vari¬ 
ous parameters by their fiducial values, we have 

ro = 0.42(a/0.1)2/®(lE/0.1)-2/®( Mbh/ 10® pc (7) 


(Eq 42, |Goodman||2003^ , where a is radiative efficiency, Ie luminosity in Eddington units, Ma 
Mach number of the viscous disc at ro, and n and Ke opacity and electron-scattering opacity, 
respectively. Hence the feeding rate to the accretion disc that eventually accretes to the SMBH 
is 

Mfeed = M,((ro/rB)'/'JB/Jc)tdi, (8) 


where the angular momentum at ro is Jo = (ro/rB)^'^^ Jb for a Keplerian disc and tdyn = Ikpc/A/Sa 
is the free-fall time at Ikpc. For our analysis, we use 


ro = O.42(MBH/lO®M0)3/25 pc, 


(9) 


bearing in mind that uncertainties are at least on the order of unity. To see how uncertainty in 
/3 affects results, we note, a 25% deviation in /3 from unity causes Mfeed in Eq ^ to change 
by a factor of 2.7, which can be compensated by adjusting each of the parameters in Eq 0 
except Mbh by a factor of 2.5 appropriately. 


3.3. Results 


We define a ratio R = 500Mfeed/SFR (SFR is the star formation rate) such that, if R is 
about unity, the observed SMBH to bulge mass ratio of ~ 0.2% (e.g., [Marconi & Hunt|||2003 


Haring & Rix 2004 would be borne out. Transformation from stellar disc(s) to a bulge is not 


addressed here. It is noted, however, that stellar discs formed from multiple gas inflows of 
inclined angles over the lifetime of a galaxy may be conducive to bulge formation. Note that 
SFR is computed directly during the simulation, whereas the SMBH accretion rate is computed 
in post-processing by evaluating Eq Figure [l] shows histories of Mfeed (blue) and R (red) 
for four random example galaxies. The most noticeable feature is that, without any intentional 
tuning, R hovers close to unity with fluctuations of order unity. 


Figureshows R as a function of redshift. We see that R increases with increasing red- 
shift from ~ 0.7 at = 0.6 - 1 to ~ 1.5 at = 3 - 4 for galaxies with Mq (green). 
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Fig. 1.— shows histories of the feeding rate Mfeed (blue) and R = 500Mfeed/SFR (red) for four 
random galaxies. The logarithm of the stellar mass for each galaxy a\. z = 0.62 is indicated at 
the top of each panel. 


with similar trends for other mass ranges. We highlight three implications. First, the observed 
SMBH to bulge ratio is readily achievable in a cosmological setting, with a slight tendency for 
more massive galaxies to have higher R. This is due to the rationing of gas supply to the 
central regions of galaxies: a small amount of gas of the lowest angular momentum feeds the 
SMBH accretion disc, while the rest builds up the stellar bulge, with the demarcation line de¬ 
termined by the accretion disc stability condition. Note that our analysis is solely based on the 
angular momentum distribution of gas that has already landed in the central Ikpc region. The 
frequency of gas inflow events into the central regions and the mass distribution of events are 
computed directly in our simulations. Second, R increases with increasing redshift, to within 
a factor of ~ 2. The trend with redshift is expected in a cosmological context, because both 
the frequency and strength of galaxy interactions increase with increasing redshift, yielding 
overall inflow gas of lower angular momentum hence a larger R at high redshift. Third, the 
smoothness of R on cosmological time scales (> lOOMyr) suggests that the dispersion of R 
is modest, around order unity, at all redshifts, consistent with the dispersion of the observed 
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Fig. 2.— shows the median of R as a function of redshift, separately for three stellar mass 
ranges M© (red), (blue) and Mq (green). The stellar mass is 

measured at the redshift in question. The vertical errorbars indicate the interquartile range, 
whereas the horizontal errorbars represent the redshift range of the bin. The red and blue 
points are horizontally slightly right-shifted for clarity of display. There are (659,2214) galaxies 
with stellar mass in the range M© for z = (3 - 4, 0.62 - 1), respectively. 

correlation locally (note that the comparison is made between computed Mfeed/SFR and ob¬ 
served Mbh/ Mbg)- Future observations at high redshift may be able to test these predictions. 
Although R is relatively smooth over cosmological time scales, the gas inflow rate varies up to 
an order of magnitude (Figure[T). The fluctuations in the inflow rate are caused by a variety of 
physical processes, including interactions between galaxies in close proximity, minor mergers 
and occasional major mergers. We have not studied in sufficient detail to ascertain whether 
secular processes play any major role. 

Is SMBFI accretion rate directly dictated by the feeding rate from galactic scales? Fig¬ 
ure]^ shows the probability distribution of feeding rate in units of Eddington rate as a function 
of Eddington ratio. The Eddington ratio is based on the assumed Mbh from the observed 
Mbh - rr relation. At 2 ; ~ 0.6 where comparisons with observations may be made, the com- 
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Fig. 3.— shows the probability distribution of feeding rate in units of Eddington rate per 
logarithmic Eddington ratio interval, as a function of Eddington ratio, in four redshift ranges, 
= 0.62 - 1 (solid red), z = I - 2 (dotted blue), = 2 - 3 (dashed green), and = 3 - 4 


(dot-dashed black) for galaxies in the stellar mass range of 10 


10.5-11 


Mq (other stellar mass 


ranges have similar properties). Also show as solid dots is the observed powerlaw distribution 
with a slope of ~ -0.6 at z ~ 0.6 from |Aird et al.| ( |2012| . The slope of the solid red curve is 
-3.3 measured for the loges range from -2.3 to -1.6 indicated by the red dashed line. 


puted distribution is steeper, computed slope -3.3 versus -0.60 observed. This indicates that 
accretion onto the SMBHs is “filtered" through physical processes operating on the accretion 
disc. This suggests that temporal correlation between AGN and star formation activities in 
individual galaxies below 30 - ISOMyr is expected to be weak, in excellent agreement with 
observations (e.g., [Hickox et al.|[2M^ . A comparison between the distribution of the feeding 
rate to the accretion disc (red curve) and that of the observed Eddington ratio (black dots) 
suggests that at ~ 0.6 accretion discs around SMBHs spend most of the time accumulating 
gas, at feeding rate below 1% Eddington ratio and that the apparent powerlaw distribution of 
Eddington ratio may be a result of superposition of AGN internal light profiles that are universal 
in shape (i.e., slope of ~ -0.6). We see that the computed feeding rate distribution shifts to 























the right ~ 0.5 dex per unit redshift, indicating that the duty cycle of luminous AGNs increases 
with redshift. 


4. Conclusions 

We have shown that, baring implausible physical fine-tuning, neither the central massive 
objects - SMBHs or NSCs - nor the stellar bulges can be regulated by blowing away the ma¬ 
jority of gas that has already landed, to explain the observed CMO-bulge relation. This leaves 
us with only one viable option. That is, the ratio of feeding rate to the nuclear region to that to 
the bulge is proportioned cosmologically. 

We test this scenario using high resolution, large-scale cosmological hydrodynamic sim¬ 
ulations without AGN feedback. Our analysis finds a proportionality, ~ 0.1 - 0.3%, between 
the feeding rate of very low angular momentum gas that can free-fall to the sub-parsec region 
to accrete to the SMBH and the star formation rate in the galaxy. There is indication that this 
ratio increases with increasing redshift to within a factor of ~ 2, suggesting that the SMBH to 
bulge ratio is nearly redshift independent, with a modest increase with redshift. We predict that 
the duty cycle of luminous AGNs increases with redshift. While SMBHs and bulges are found 
to coevolve on > 30 - 150Myr time scales, there is indication that, on smaller time scales, 
the SMBH accretion and star formation may be less or not correlated, which is likely due to 
variations of AGN activities on smaller time scales dictated by physics of accretion disc. 

While our analysis disfavor internal regulation in terms of blowing gas away with the re¬ 
quired proportionality, “random" internal regulation by blowing some gas away without the said 
proportionality is not ruled out and in fact may be common, manifested as galactic superwinds 
or AGN winds. Nor do we disfavor feedback processes that control the overall amount of cold 
gas supply, termed "global feedback". Global feedback reflects the collective effects of stel¬ 
lar evolution (supernovae, winds, etc) and SMBH accretion (winds, radio jets, etc) as well as 
gravitational shock heating due to structure formation and photoionization heating, among oth¬ 
ers. They impact the thermodynamical state of the interstellar, circumgalactic and intergalactic 
medium. We emphasize that, even if global feedback controls the overall cold gas supply and 
its temporal distribution on cosmological time scales, it is not responsible for the proportional 
growth of SMBHs and galaxies. 

An implication is that the distinction between forming a NSC or SMBH may hinge on the 
existence of a massive enough initial black hole seed. Thus, the demarcation bulge mass of 
Mbgo = (3-4) X10^° Mq is suggestive that only the progenitors of the massive enough galax¬ 
ies have formed massive black hole seeds at some high redshift, with less massive galaxies 
seeded by NSCs or neither. Subsequently, those with initial massive black hole seeds are able 
to accrete the infallen gas and grow to SMBHs over time, whereas those without massive black 
hole seeds turn the infallen gas in the nuclear regions into stars to grow the NSCs. Let us sup¬ 
pose that CMOS of initial mass Mcmomu created at some high redshift in dwarf galaxies have 
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migrated to the centers of larger galaxies to serve as central seeds. The rationed gas supply 
would then yield final Mcmo/Mbg = (a MBG + McMo,init)/( MBG + ^cMo,init)- Thus, for those 
galaxies lacking significant, subsequent growth of the CMO, i.e., a Mbg is not much greater 
than McMo,init, the CMO-bulge mass scaling relation will be sublinear, which may explain the 
observed shallower scaling relation between NSCs and bulges at the low end of bulge mass 
(e.g., [Erwin & Gadotti] 201^ [Leigh et aL]|2012| [Scott & Graham|[20T3| [den Brok et al.|2014[ ). 
Galaxies with a massive initial black hole seed may form a NSC as well, consistent with obser¬ 
vations (e.g., [Seth et al.[[2008| [Gonzalez Delgado et al.S2008^ , although the stellar component 
in the vicinity of an SMBH may be altered by subsequent, additional processes, such inspiral 
of another SMBH (e.g., [Milosavljevic et al.[|2002^ . 

This study is related to [Escala| ( |2006[ [2007| , who studied gas accretion processes sur¬ 
rounding the SMBH; we explicitly avoid detailed accretion physics by focusing on the amount 
of mass that enters the “feeding" zone to the SMBH. This work reaches conclusions similar to 
that of [Angles-Alcazar et'ar]|2015| with respect to the Mbh - Mbg ratio, with a contrasting 
difference on the role of feedback. While [Angles-Alcazar et aL] ( |2015) requires that only a 
small fraction of the gas at subparsec scales is actually accreted by the SMBH, with the rest 
lost to winds and outflows, we suggest that the gas disc beyond the Toomre unstable radius 
is instead consumed by star formation, without requiring blowing away most of the gas by the 
SMBH. 


I am indebted to an anonymous referee for the most detailed, cogent, critical yet civilized 
reports, which have immensely helped improve the presentation and clarify numerous issues. 
I thank Dr. Guangtun Zhu for very helpful discussion. This work is supported in part by grant 
NASA NNX11AI23G. 
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